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Foreword 


One  of  the  Major  problems  in  using  satellite  data  to 
study  mesoscale  ocean  features  is  that  of  atmospheric 
backscattering  of  sunlight,  which  has  a  masking  effect 
on  the  data. 

This  report  deals  with  a  new  atmospheric  correction 
method,  used  with  the  Nimbus-7  Coastal  Zone  Color 
Scanner  (CZCS),  that  effectively  removes  atmospheric 
radiance  from  the  satellite  images.  The  method  em¬ 
ployed  is  Principal  Components  Analysis,  a  statistical 
tool  used  in  the  field  of  pattern  recognition. 

This  method  is  easy  to  use  and  allows  CZCS  data  to  be 
incorporated  into  routine  naval  environment-oriented 
operations.  This  report  presents  four  examples  for 
three  ocean  areas. 


J3. 

G.  T.  Phelps,  Captain,  USN 
Commanding  Officer,  NORDA 


Executive  Summary 


The  Nimbus-7  Coastal  Zone  Color  Scanner  (CZCS)  has  the 
unique  potential  to  remotely  sense  mesoscale  ocean 
features  through  warm,  humid  atmospheres  that  are 
opaque  to  thermal  infrared  sensors.  The  major  obstacle 
to  the  use  of  these  data  is  the  masking  effect  of  sun¬ 
light  backscattered  by  the  atmosphere.  This  study 
details  a  new  atmospheric  correction  method,  employing 
principal  component  analysis  techniques,  that  effec¬ 
tively  removes  the  atmospheric  radiance  from  CZCS 
images  of  the  ocean.  The  method  has  the  advantage  of 
ease  of  use  and,  hence,  allows  the  CZCS  data  to  be 
incorporated  into  routine  naval  environmentally- 
oriented  operations.  The  study  shows  the  method  to  be 
especially  effective  in  the  descriptive  analyses  of 
mesoscale  oceanic  phenomena.  Four  examples  of  the  use 
of  the  technique  are  presented  for  three  different 
ocean  areas . 
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Scanner  Data  to  Define  Mesoscale  Ocean  Features  through  a  Warm,  Humid  Atmosphere 


1.  Introduction 

Oceanographers  utilizing  satellite  data 
to  study  ocean  mesoscale  processes  fre¬ 
quently  use  infrared  (IR)  sensors  that 
operate  in  the  8,000—13,000  nanometer 
(nm)  range  to  produce  imagery  of  the 
radiation  temperatures  of  the  ocean  sur¬ 
face  (Fig  •  1).  Unfortunately,  these 

thermal  images  do  not  represent  exact 
sea  surface  temperatures  because  the 
data  include,  in  addition  to  ocean  radi¬ 
ation,  the  absorption  and  radiation 
effects  of  the  intervening  atmosphere* 
The  atmospheric  constituent  that  con¬ 
tributes  the  most  to  corrupting  ocean 
radiation  at  infrared  wavelengths  is 
water  vapor*  For  a  cool,  dry  atmosphere, 
the  effects  due  to  water  vapor  are  suf¬ 
ficiently  small  so  that  sea  surface  tem¬ 
perature  patterns  are  readily  observed 
in  IR  imagery.  In  these  cases,  descrip¬ 
tive  analyses  may  be  performed  without 
the  need  for  water  vapor  correction* 
However,  as  the  amount  and  temperature 
of  the  atmospheric  water  vapor  increa¬ 
ses,  the  ocean  thermal  patterns  in  the 
imagery  become  increasingly  obscure. 
Figure  2,  a  summer  image  of  the  Gulf  of 
Mexico,  is  a  good  example  of  the  ocean 
viewed  by  an  IR  sensor  through  a  humid 
atmosphere.  This  figure  may  be  compared 
with  the  winter  IR  image  of  the  Gulf  in 
Figure  1  to  show  the  lack  of  observable 
ocean  thermal  expression  during  the  sum¬ 
mer  months. 

A  graphic  demonstration  of  the  effects 
of  the  atmosphere  in  infrared  imagery  is 
shown  in  Figure  3.  This  NOAA-7  Advanced 
Very  High  Resolution  Radiometer  (AVHRR) 
image  shows  a  strong  atmospheric  water 
vapor  front  present  over  the  Gulf  of 
Mexico  on  19  June  1979.  The  image  shows 


more  oceanographic  details  on  the  dry 
air  side  of  the  front  (i.e.,  in  the 
northeastern  Gulf)  than  in  the  more 
humid  southwestern  side. 

The  interference  of  the  humid  atmos¬ 
phere,  such  as  that  shown  in  these 
figures,  generally  renders  satellite 
infrared  sensing  of  the  ocean  useless 
over  the  tropical  oceans  during  all 
seasons,  and  of  limited  value  at  middle 
latitudes  over  significant  portions  of 
the  year. 

Visible  and  near-infrared  portions  of 
the  electromagnetic  spectrum  are  not  as 
severely  contaminated  by  atmospheric 
water  vapor  effects  as  are  infrared 
wavelengths  (aerosols  are  the  major 
problem  at  visible  wavelengths).  There¬ 
fore,  multi-channel  visible  and  near- 
infrared  spectrometers  such  as  the 
Nimbus-7  Coastal  Zone  Color  Scanner 
(CZCS)  offer  the  possibility  of  imaging 
ocean  features  through  warm,  humid 
atmospheres  and  may,  therefore,  open  the 
tropical  and  summer  mid-latitude  oceans 
to  studies  using  satellite  remote  sens¬ 
ing  techniques. 

The  Nimbus-7  CZCS  is  a  six-channel  scan¬ 
ning  radiometer.  These  channels  measure 
the  earth’s  radiance  in  spectral  bands 
centered  at  443  nm,  520  nm,  550  nm,  670 
nm ,  750  nm ,  and  11,500  nm .  The  first 
four  channels  are  in  the  visible,  the 
fifth  is  in  the  near  (reflective)  infra¬ 
red,  and  the  sixth  is  in  the  thermal 
(emissive)  infrared  portions  of  the 
electromagnetic  spectrum.  Further  de¬ 
tails  on  the  sensor  may  be  obtained  from 
Hovis  et  al.  (1980). 

There  are  three  major  questions  that 
must  be  answered  to  determine  if  this 


Figure  1.  NOAA-7  thermal-IR  (11.0/i)  image  of  the  Gulf  of  Mexico  for  10  December 
1981.  The  image  clearly  shows  the  warm  (dark  shades)  Loop  Current  flowing  into 
the  Gulf  through  the  Yucatan  Straits  and  the  seasonally  cooled  light  shades  shelf 
waters  adjacent  to  the  landmasses,  especially  Florida.  This  image  is  typical  of  IR 
remote  sensing  through  a  cool,  dry  atmosphere. 


potential  of  Nimbus-7  CZCS  imagery  can 
be  realized  for  naval  operational  use: 

1)  Can  ocean  color  features  be  imaged 
from  space  using  the  CZCS  through  the 
backscattered  radiance  of  the  earth’s 
atmosphere? 

2)  Are  the  ocean  color  features  seen  in 
the  CZCS  imagery  associated  with  meso- 
scale  oceanic  structures  such  as  fronts 
and  eddies? 

3)  If  ocean  color  features  can  be  imaged 
using  CZCS  data,  and  they  do  represent 
ocean  structure,  can  a  routine  imaging 
method  be  derived  that  is  suitable  for 
Navy  use? 


An  example  for  the  rationale  for  Ques¬ 
tion  1  is  illustrated  in  Figure  4.  This 
raw  CZCS  image  of  the  Grand  Banks  region 
of  the  North  Atlantic  does  not  show  any 
major  surface  structure  in  this  oceano¬ 
graphically  dynamical  area.  Any  varia¬ 
bility  in  water  color  in  the  data  is 
evidently  so  subtle  as  to  be  completely 
masked  by  the  atmospheric  path  radiance. 
Atmospheric  path  radiance  refers  to  the 
sunlight  reflected  to  the  satellite  from 
constituents  in  the  air  between  the  sat¬ 
ellite  and  the  ocean.  Normally  more  than 
90%  of  the  radiance  observed  from  satel¬ 
lite  altitudes  over  the  ocean  is  due  to 
atmospheric  path  radiance.  Thus,  to 
achieve  reasonable  measurements  of  ocean 
radiance,  it  is  necessary  to  remove  this 
effect. 
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Figure  2.  TIROS-N  thermal-IR  (11.0  /i )  image  of  the  Gulf  of  Mexico  for  24  June 
1  979.  Note  that  there  is  no  discernible  evidence  of  the  Loop  Current  or  other 
oceanographic  features.  This  image  is  typical  of  IR  remote  sensing  through  a  warm, 
humid  atmosphere. 


The  most  widely  accepted  method  for  cor¬ 
recting  for  atmospheric  path  radiance  in 
CZCS  data  utilizes  a  model  based  on 
single-scattering  theory  (Gordon,  1978). 
The  single-scattering  model  requires 

•  several  assumptions  about  both  the  scat¬ 
tering  and  absorption  properties  of  the 
atmosphere  and  ocean.  Although  these 
assumptions  are  valid  in  many  cases, 
their  use  in  correcting  CZCS  data  can 
lead  to  significant  uncertainties  in 

•  other  cases. 


This  report  will  demonstrate  some  of  the 
problems  and  uncertainties  encountered 
in  applying  the  Gordon  atmospheric  cor¬ 
rection  and  will  present  an  alternative 
correction  method.  The  method — Principal 


Components  Analysis  (PCA) ,  a  statistical 
tool  widely  used  in  the  field  of  pattern 
recognition — avoids  many  of  the  problems 
associated  with  the  single— scattering 
model  approach.  However,  it  should  be 
noted  that  at  the  stage  of  development 
presented  here,  the  PCA  method  does  not 
result  in  radiance  values  but  generates 
an  atmospherically  corrected  image  in  an 
abstract  mathematical  space.  For  this 
reason  the  PCA  method  is  not  useful 
where  absolute  spectral  radiance  values 
are  required.  Thus,  it  is  not  presented 
as  a  replacement  for  the  single¬ 
scattering  model  methodology;  rather, 
PCA  is  useful  in  operationally  producing 
CZCS  imagery  of  ocean  features  of  a 
quality  that  allows  full  descriptive 
interpretation. 
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Figure  3.  TIROS-N  thermal-IR  (11.0  pt )  imagery  of  the  Gulf  of  Mexico  on  19  June 
1979  showing  an  atmospheric  temperature/humidity  front.  Note  the  "clean"  appear¬ 
ance  of  the  northeastern  Gulf  compared  with  the  "foggy"  appearance  of  the  south¬ 
west  quadrant  of  the  figure.  The  boundary  between  the  two  atmospheric  regimes  is  ® 

quite  sharp  and  splits  the  Loop  Current  in  half. 
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Figure  4.  Color  image  formed  from  CZCS  channels  1,  2,  and  3  of  orbit  2811  on 
15  May  1979  at  1456  Z.  Image  is  from  raw  data  with  no  enhancements  applied. 

•  Box  outlines  the  area  used  for  PCA. 


2.  Atmospheric  Correction  Techniques: 

The  Single-Scattering  Method 

•  The  following  overview  of  the  CZCS 

atmospheric  correction  problem  will  help 
identify  some  of  the  sources  of  uncer¬ 
tainty  in  the  single-scattering  model 
methodology . 


t(VLw(V =  MV  ■  MV  "  S(W 
{L^V  "  W  -tW(V}  (2) 


where  S(A  ,A  )  =  L  (A  )/L  (A  )  . 

l  Z  a  I  a  Z 


•  2. 1  The  Single-Scattering 

Approximation 

Single-scattering  theory,  as  the  name 
implies,  ignores  multiple  scattering. 
The  theory  assumes  that  a  photon  travel- 

•  ing  the  atmospheric  path  from  space  to 
the  earth’s  surface  is  unperturbed  or 
experiences,  at  most,  one  interaction 
with  an  atmospheric  constituent.  This 
theory  is  reasonable  for  optically  thin 
atmospheres  but  begins  to  break  down  for 

•  atmospheres  that  are  more  optically 
dense.  The  advantage  of  single¬ 
scattering  theory  is  that  the  Rayleigh 
(molecular)  and  aerosol  (particulate) 
radiance  contributions  to  the  backscat- 
tered  radiance  of  the  atmosphere  are 

•  additive.  Thus,  the  observed  radiance  of 
the  ocean  (and  atmosphere)  as  viewed 
from  space  is  given  by 

L(A)  =  L  (A)  +  L  (A)  +  t (A)L  (A)  (1) 

r  a  w 

where  L(A)  ,  L  (A),  L  (A),  and  L  (A)  are 
r  a  w 

the  observed,  Rayleigh,  aerosol,  and 
ocean  radiances  respectively  at  wave¬ 
length  X,  and  t(A)  is  the  transmittance 

•  of  the  atmosphere.  In  Eq.  (1),  L(A)  is 
observed  by  the  CZCS  radiometer  and 
L^(A)  can  be  calculated  with  reasonable 

accuracy  (Frohlich  and  Shaw,  1980).  How¬ 
ever,  L  (A)  and  L  (A)  are  inseparable, 
w  a  y 

0  so  a  direct  determination  of  t(  A )L  (X) 

is  not  possible.  w 

To  solve  for  L  (A)  ,  Gordon  (1978)  uses 

two  wavelengths  simultaneously.  Taking 
two  equations  like  Eq.  (1),  one  at  \2, 

•  and  the  other  at  Xj,  the  solution  of 
simultaneous  equations  yields 


Gordon  (1981)  gives  the  following 
expression  for  S(X1>  X2)* 

s<W  *  C<W  [VV'VV] 

exp  {-  iw-vyj 

(3) 


where  Tq(X)  is  the  ozone  optical  thick¬ 
ness,  F  ( X )  is  the  extraterrestrial 

solar  irradiance,  n  and  u  are  the 

o  s 

cosines  of  the  viewing  angle  and  solar 
zenith  angle  respectively,  and  e(X]L,  X2) 

is  related  to  the  aerosol  optical  thick¬ 
ness  x^(X)  and  the  single-scattering 

albedo  wq  through 

e<W  =  VVVW^WV* (4) 

To  solve  Eq.  (2)  for  t(A, )L  (A.),  one 

1  w  1 

must  make  assumptions  concerning  L  (\n) 

w  2 

and  e(X1,X2). 

2.2  Assumptions  About  Oceanic 
Spectral  Radiance 

The  term  in  Eq.  (2)  for  water  radiance 
at  X2>  hanc^-ed  by  selecting 

X2  such  that  Lw(X2)  =  0;  thus,  this  tern 

can  be  ignored.  For  CZCS  data  the  red 
channel  (670  nm)  is  used  as  ^ .  This 

channel  is  chosen  because  the  absorption 
coefficient  of  sea  water  is  higher  at 
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this  wavelength  than  at  wavelengths 
sensed  by  the  other  CZCS  visible  chan¬ 
nels.  Thus,  relatively  small  backscatter 
from  the  ocean  can  be  expected.  Figure  5 
shows  upwelled  spectral  radiance  data 
measured  just  above  the  ocean  surface  at 
several  stations  in  the  Gulf  of  Mexico 
during  December  1978.  The  data  shows 
that  the  amount  of  upwelled  radiance 
does  fall  off  drastically  at  red  wave¬ 
lengths  (actually,  most  of  the  observed 
radiance  at  670  nm  seen  in  the  Figure  5 
data  is  skylight  reflecting  off  the  sea 
surface  rather  than  subsurface  backscat¬ 
ter)  .  Thus,  the  assumption  that  ocean 
radiance  in  CZCS  channel  4  is  negligible 
seems  to  be  reasonable  for  the  deep 
ocean  (i.e.,  low  turbidity)  cases.  How¬ 
ever,  this  assumption  is  not  valid  for 
coastal  regions  where  significant  back¬ 
scatter  does  occur  at  670  nm.  For  exam¬ 
ple,  Neville  et  al .  (1980),  reporting  on 


Figure  5.  Spectral  radiance  of  the  ocean 
measured  from  just  above  the  sea  sur¬ 
face  at  several  locations  in  the  Gulf  of 
Mexico  during  November  and  December 
1978. 


a  study  of  CZCS  data  for  Lake  Ontario, 
cite  appreciable  reflectance  at  670  nm 
as  the  reason  for  the  poor  performance 
of  the  single-scatter  method  in  their 
study.  To  help  offset  this  problem,  an 
iterative  algorithm  has  been  developed 
by  Smith  and  Wilson  (1980). 

2.3  Assumptions  About  Uncertainty 
in  Spectral  Distribution  of 
Aerosol  Backscatter 


The  factor  e(X^,  \^)  in  Eq.  (5)  repre- 

sents  the  spectral  variability  in  the 
aerosol  phase  function.  However,  infor¬ 
mation  on  atmospheric  aerosol  phase 
functions  is  sparse.  For  simplicity  it 
can  be  assumed  that  the  aerosol  optical 

thickness  T  (\)  ,  varies  with  wavelength 
a 

according  to  the  power  law 


t  a)  =  t  (X  )  (A  /A)' 

a  a  o  o 


(5) 


where  n  is  referred  to  as  the  Angstrom 
exponent  (Angstrom,  1961). 


To  make  an  atmospheric  correction  to 
CZCS  data,  one  must  either  measure  n  at 
the  time  of  satellite  overpass  or  assume 
a  value  for  n.  Gordon  (1978)  assumes 
that  the  aerosol  phase  function  depends 
weakly  on  wavelength  over  the  visible 
part  of  the  spectrum,  i.e.,  that  the 
value  of  n  is  approximately  zero.  How¬ 
ever,  this  is  not  strictly  true  in  all 
cases  and  the  Angstrom  exponent  for 
atmospheric  aerosols  has  been  found  to 
range  between  0  and  2  (Robinson,  1966). 
Tomasi  and  Prodi  (1982)  in  measuring 
atmospheric  turbidity  in  the  Red  Sea, 
Indian  Ocean,  and  Somalian  Coast  found  n 
values  that  ranged  between  0.2  and  0.9. 


Measurements  of  the  aerosol  optical 
depth  of  the  atmosphere  made  in  November 
and  December  1978  in  the  Gulf  of  Mexico 
are  shown  in  Figure  6,  along  with  power 
law  fits  of  the  measured  values.  The 
figure  shows  that  the  power  law  fit  to 
aerosol  optical  depth  is  only  approxi¬ 
mate  and  that  the  Angstrom  exponent 
actually  varied  from  day  to  day;  ranging 
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Figure  6.  Aerosol  optical  depth  of  a  ver¬ 
tical  atmospheric  path  measured  in  the 
Gulf  of  Mexico  onboard  the  USNS  LYNCH 
on  three  different  days  in  December 
1978.  Dotted  lines  are  power  law  fits  to 
the  measured  data. 

from  0  to  1.9  with  a  mean  value  of  1.2. 
This  mean  is  close  to  the  1.3  value 
often  used  for  continental  regions 
(Robinson,  1966).  Because  of  this  type 
of  aerosol  variability,  large  errors  in 
correction  terras  can  result  when  670  nn 
values  are  extrapolated  into  the  green 
and  blue  regions  based  on  some  assumed 
value  for  n. 

Gordon  and  Clark  (1981)  propose  using 
spectral  radiance  values  from  a  clear 
water  portion  of  a  CZCS  image  to  derive 
values  for  the  Angstrom  exponent.  They 
state  that  good  results  are  possible  if 
a  water  area  with  less  than  0.25  mg/m3 
of  chlorophyll-a  can  be  found  within  a 
given  scene. 

The  single-scatter  modeling  approach  to 
atmospheric  correction  may,  therefore. 


be  summarized  as  follows:  Each  of  the 
four  CZCS  channels  are  corrected  for  the 
Rayleigh  backscatter  of  the  atmosphere. 
Since  the  ocean  is  assumed  to  be  "black" 
at  670  nm,  any  radiance  remaining  in 
channel  4  after  the  Rayleigh  correction 
is  attributed  to  atmospheric  aerosols. 
The  Rayleigh-corrected  red  (670  nm) 
channel  is,  therefore,  the  aerosol  cor¬ 
rection  term  that  is  applied  to  the 
remaining  three  channels.  However,  be¬ 
fore  application,  the  correction  term  is 
differentially  weighted  for  each  of  the 
three  channels.  The  weighting  factors 
are  somewhat  subjective  and  require 
knowledge  of  the  solar  spectral  irradi- 
ance,  the  ozone  optical  depth,  and  the 
Angstrom  exponent  for  the  atmospheric 
aerosols  for  each  case. 

Thus,  the  corrected  CZCS  image  may  dif¬ 
fer  significantly  from  the  actual  ocean 
spectral  radiance  for  three  important 
reasons : 

•  the  single-scattering  approximation 
may  be  poor  because  the  atmosphere  is 
optically  thick; 

•  the  ocean  is  not  really  "black"  at  670 
nm; 

•  an  incorrect  value  has  been  used  for 
the  Angstrom  exponent. 

2.4  The  Single-Scattering  Model 
Applied  to  a  Test  Case 

A  simplified  form  of  the  single¬ 
scattering  model  applied  to  an  ocean 
data  set  will  both  illustrate  the  essen¬ 
tial  elements  of  this  established  method 
and  serve  as  a  reference  for  evaluating 
the  proposed  statistical  method.  The 
data  used  is  from  the  Grand  Banks  CZCS 
image  shown  in  Figure  4.  The  CZCS  ther¬ 
mal  infrared  image  (channel  6)  corre¬ 
sponding  to  the  CZCS  visible  data  is 
shown  in  Figure  7.  This  area,  one  of  the 
most  dynamic  regions  in  the  oceans, 
forms  the  confluence  of  three  major 
current  systems:  the  Gulf  Stream,  the 
Labrador  Current,  and  the  North  Atlantic 
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Figure  7.  Image  of  CZCS  channel  6  (infrared)  data  coincident  with  Figure  4.  Labels 
on  image  identify  major  oceanographic  regimes. 


Current.  Note  that  although  the  surface 
thermal  gradients  created  by  the  encoun¬ 
ter  of  the  thermally  different  currents 
are  well-delineated  in  the  IR  imagery, 
no  corresponding  features  are  discerni¬ 
ble  in  the  visible  color  image.  Either 
there  is  no  color  signature  associated 
with  these  divergent  current  systems,  or 
the  ocean  surface  upwelling  radiance  is 
obscured  by  the  more  dominant  atmospher¬ 
ic  backscatter.  The  data,  therefore, 
present  an  ideal  case  to  demonstrate  the 
two  atmospheric  correction  procedures. 

The  single-scattering  model  atmospheric 
correction  was  performed  on  the  data 
following  the  general  methodology  previ¬ 
ously  described.  Clear  water  radiances 


(Gordon  and  Clark,  1981)  were  used  to 
derive  an  Angstrom  exponent  of  2.8  for 
this  case.  Solar  irradiance  values,  FQ 

(X),  were  assumed  to  be  184.63,  183.57, 
185.01,  and  153.13  in  CZCS  channels  1, 
2,  3,  and  4,  respectively  (Austin, 

1981).  The  values  for  ozone  transmit¬ 
tance  of  the  atmosphere,  j  (X),  were 

taken  to  be  0.0009,  0.0146,  0.0266,  and 
0.0138  from  the  L0WTRAN-5A  code  (Kneizys 
et  al.,  1980)  for  a  mid-latitude  summer 
atmosphere.  Solar  and  satellite  zenith 
angles  were  held  constant  at  their 
values  in  the  box  drawn  in  Figure  4. 
Solar  and  satellite  zenith  angles  at 
that  location  are  23.5°  and  25°.  Not 
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Figure  8.  Color  image  formed  from  atmospherically  corrected  CZCS  channels  1,  2, 
and  3.  Atmospheric  correction  was  by  means  of  a  simplified  single-scattering 
atmospheric  model. 


permitting  these  angles  to  vary  spatial¬ 
ly  results  in  Rayleigh  correction  errors 
of  15%  at  the  left  edge  of  the  image  in 
the  figure  and  25%  at  the  right  edge. 
(The  limb  effects  seen  in  the  corrected 
data  would  largely  disappear  if  angular 
dependence  was  included  in  the  Rayleigh 
processing . ) 

When  the  single- scattering  correction  is 
applied,  the  corrected  radiance  values 
in  channel  1  (443  nm)  become  negative. 
Since  this  is  physically  impossible,  it 
implies  that  the  correction  method,  or 
at  least  the  values  of  the  method  param¬ 
eters,  is  erroneous.  Negative  corrected 
radiances  could  result  from  several 
sources  including  uncertainties  in  solar 
spectral  irradiance  or  sensor  calibra¬ 
tion.  In  this  case  the  Angstrom  exponent 
value  of  2.8,  derived  from  clear  water 
radiances,  seems  suspiciously  large 
(since  the  Angstrom  exponent  is  general¬ 
ly  considered  to  range  from  0  to  2). 
There  is  no  obvious  explanation  why  the 
clear  water  radiances  gave  such  a  large 
value,  nor  is  there  ground  truth  or 
other  collaborating  data.  Prelaunch  cal¬ 
ibration  data  are  used  to  convert  raw 
data  to  radiance  values.  The  prelaunch 
and  in-flight  calibration  to  GZCS  is 
presently  a  matter  of  considerable 
uncertainty.  It  could  be  that  poor  cali¬ 
bration  data  led  to  the  large  n  value. 
If  the  clear  water  radiance  approach  is 
abandoned  and  a  mid-range  Angstrom  expo¬ 
nent  of  1  is  assumed,  the  following 
values  S(A^,  A^)  result: 

•  S(443 ,670)  =  1.80; 

•  S(520, 670)  -  1.56; 

•  S(550 ,670)  =  1.44. 

The  atmospheric  correction  of  the  CZCS 
image  is  repeated  using  these  values. 
The  channel  1  radiances  now  come  out 
positive.  A  subjective  visual  examina¬ 
tion  of  the  corrected  imagery  indicates 
that  a  reasonable  correction  of  the  data 
has  been  made. 

Figure  8  is  a  color  composite  of  the 
three  channels  of  atmospherically  cor¬ 
rected  data  that  shows  an  ocean  differ¬ 
ent  than  that  shown  in  the  CZCS  image  of 


Figure  4.  The  correction  appears  to  have 
effectively  stripped  away  the  atmospher¬ 
ic  path  radiance  revealing  oceanic  color 
patterns  similar  to  the  general  features 
in  the  thermal  image  of  Figure  7  . 

3.  Atmospheric  Correction  Techniques: 
Principal  Components  Analysis  Method 

To  overcome  some  of  the  single-scatter 
method  problems  in  the  case  just  de¬ 
scribed,  the  atmospheric  correction  of 
CZCS  imagery  can  be  approached  as  a 
statistics  problem.  With  this  approach, 
assumptions  about  physical  properties 
that  are  part  of  the  Gordon  (1978)  meth¬ 
od,  such  as  that  the  ocean  is  black  at 
670  nm,  or  that  the  aerosol  phase  func¬ 
tions  vary  only  weakly  with  wavelength, 
are  replaced  with  assumptions  about  the 
statistical  properties  of  the  data.  With 
this  method  an  assumption  is  made  that 
there  is  an  area  of  the  image  where  most 
of  the  variance  in  the  multi-spectral 
measurement  space  is  attributable  to 
aerosol  variability.  Two  important  cri¬ 
teria  of  the  area  where  this  statistical 
property  exists  are:  1)  the  area  should 
have  a  high  variability  in  aerosol  con¬ 
centration,  and  2)  the  ocean  be  nearly 
uniform  in  color.  For  this  reason  coast¬ 
al  waters,  with  their  typically  high 
water  color  variability,  are  avoided. 
Likewise,  areas  with  clouds  are  avoided 
since  the  brightness  of  a  cloud  would 
saturate  the  eight-bit  dynamic  range  of 
the  CZCS.  Such  saturation  introduces  a 
nonlinearity  into  what  is  basically  a 
linear  transformation  technique. 

3.1  Principal  Components  Analysis 

Principal  Components  Analysis  (PCA)  is 
most  frequently  used  in  pattern  recogni¬ 
tion  problems  to  reduce  the  dimensional¬ 
ity  of  a  measurement  space.  For  example, 
it  is  possible  to  transform  10-channel 
raultispectral  imagery  into  only  three 
channels,  and  yet  preserve  95%  or  more 
of  the  "information"  present  in  the 
original  10-channel  data.  In  PCA  "infor¬ 
mation"  is  synonymous  with  variance; 
i.e.,  PCA  is  a  dimensionality  reduction 
tool  that  preserves  the  maximum  possible 
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variance.  In  some  situations,  the  vari¬ 
ance  present  may  be  due  to  noise  rather 
than  information,  so  an  analyst  should 
know  the  statistical  properties  of  the 
data  before  choosing  to  use  PCA.  In  the 
present  case,  the  image  subarea  selected 
for  PCA  has  been  chosen  to  have  the 
assumed  statistical  properties.  If  these 
desired  properties  do  exist  in  the  data, 
then  the  same  characteristics  of  PCA 
that  make  it  useful  for  dimensionality 
reduction  make  it  useful  for  atmospheric 
corrections . 

From  a  mathematically  nonrigorous  point 
of  view,  PCA  is  the  linear  transforma¬ 
tion  of  a  given  set  of  variables  to  a 
new  set  of  variables  called  principal 
components.  This  new  set  of  variables  is 
arranged  in  a  sequence  such  that  each 
succeeding  component  is  that  linear  com¬ 
bination  of  the  original  variables  which 
has  the  maximum  variance  of  all  possible 
combinations  under  the  constraint  that 
each  component  be  uncorrelated  with,  and 
orthogonal  to,  all  preceeding  compo¬ 
nents.  (More  details  may  be  found  in 
Anderson,  1958,  or  Cooley  and  Lohnes, 
1971.)  It  can  be  shown  mathematically 
that  the  principal  components  defined  in 
this  manner  are,  in  fact,  the  eigen¬ 
vectors  of  the  covariance  matrix. 
Furthermore,  it  can  be  shown  that  the 
variance  of  the  i-th  principle  component 
is  identically  equal  to  the  i-th  eigen¬ 
value  of  the  covariance  matrix. 

The  variance-ordering  properties  of  the 
principal  component  sequences  and  the 
fact  that  most  of  the  variance  in  the 
selected  subarea  of  the  image  results 
from  aerosol  variability  are  the  basis 
for  the  removal  of  the  aerosol  effects. 
If  the  first  principal  component  maxi¬ 
mizes  variance,  then  this  component  cap¬ 
tures  the  aerosol  contribution  to  the 
total  variance.  In  other  words,  because 
the  selected  subarea  is  bland  oceano¬ 
graphically,  the  orientation  of  the 
first  principal  component  will  be  paral¬ 
lel  to  the  aerosol  variability.  The 
first  principal  component  can,  there¬ 
fore,  be  discarded  to  make  the  aerosol 
correction.  In  the  CZCS  image  used  to 


illustrate  the  technique,  the  original 
measurement  space  is  four-dimensional 
(CZCS  channels  one  through  four).  PCA 
results  in  a  transformation  to  a  new 
four-dimensional  space.  If  the  first  new 
dimension  (or  first  component)  is 
removed  (i.e.,  projected  onto  the  other 
three  in  mathematical  terms),  the  result 
is  a  three-channel  image,  with  the 
obscuring  effect  of  atmospheric  aerosols 
removed.  Because  of  this  removal,  the 
resulting  three-channel  image  will  show 
the  oceanic  portion  of  the  variance  that 
was  hidden  in  the  original  image. 

3.2  Statistics  of  the  Test  Case  Data 


Statistics  for  the  selected  subarea  in 
the  Grand  Banks  data  set  in  Figure  4 
(boxed  in  the  image)  are  given  (p.16), 
along  with  the  principal  components  and 
the  eigenvalues  of  the  covariance  ma¬ 
trix.  As  previously  mentioned,  the  prin¬ 
cipal  components  are  formed  from  linear 
combinations  of  the  original  spectral 

bands.  If  x,,  .  is  the  i-th  subarea  pixel 
ij 

from  channel  j,  then  the  k-th  principal 
component,  P  can  be  written 


(6) 


The  coefficients  a 


kj* 


therefore,  define 


the  transformation  from  spectral  meas¬ 
urement  space  to  principal  component 
space.  The  statistics  below  are  given  in 
vector  form  where  the  vector  components 
represent  each  of  the  four  spectral 
bands.  This  data  has  not  been  radiomet- 
rically  calibrated,  so  the  values  given 
are  in  units  of  raw  digital  counts. 


Note  that  the  eigenvalues  indicate  that 
the  first  principal  component  (the 
atmospheric  aerosol  component)  in  this 
case  account  for  96.6%  of  the  total 
variance  within  the  box.  This  means  that 
only  3.4%  of  the  variance  is  of  oceanic 
origin,  which  explains  why  oceanic  fea¬ 
tures  could  not  be  seen  in  the  raw  data 
in  Figure  4. 


Figure  9.  Comparison  of  the  aerosol  corrections  for  the  two  methods.  The 
upper  figure  is  the  aerosol  correction  from  single-scattering  theory,  and  the 
lower  figure  is  the  aerosol  correction  from  the  PCA  method. 
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Case  I  Grand  Banks 
PCA  Box  Statistics 


mean  =  162,  149,  151,  148 

std  dev  =  6.57,  8.64,  10.13,  17.08 

eigenvalues  =  0.966,  0.028,  0.004,  0.001 


kj 


0.265  0.381  0.451  0.761 
-0.731  -0.382  -0.159  0.541 
-0.529  0.199  0.743  -0.356 
0.337  -0.817  0.466  0.015 


The  a^_.  values,  particularly  the  weights 

for  the  third  principal  component,  are 
also  interesting.  The  coefficients 
-0.529,  0.199,  and  0.743  for  channels  1, 
2,  and  3,  respectively  (row  3  of  the 
matrix) ,  indicate  that  this  component 
shows  changes  associated  with  chloro¬ 
phyll  concentrations.  This  statement  is 
based  on  the  fact  that  the  component  is 
composed  of  large  contributions  from 
channels  1  and  3  that  are  opposite  in 
sign  plus  a  relatively  small  contribu¬ 
tion  from  channel  2.  This  weighting 
exhibits  the  well-known  "hinge  point" 
signature  of  chlorophyll-a  (Duntley  et 
al. ,  1974).  In  a  deep  ocean  location 

like  the  Grand  Banks,  chlorophyll  is 
normally  the  primary  contributor  to 
ocean  color,  so  it  would  be  expected 
that  this  third  component  will  probably 
show  more  oceanic  structure  than  will 
the  other  components. 


correction  factor  using  the  single¬ 
scattering  method).  Visual  comparison  of 
the  two  indicates  that  an  aerosol  cor¬ 
rection  based  on  the  first  principal 
component  will  be  very  similar  to  a  cor¬ 
rection  based  on  Rayleigh  corrected 
channel  4. 

As  expected,  the  third  principal  compo¬ 
nent  shows  the  most  features  of  apparent 
oceanographic  origin.  Figure  10  is  an 
image  of  the  third  principal  component, 
which  accounted  for  0.4%  of  the  total 
variance.  Comparison  of  Figure  10  with 
the  IR  data  (Fig.  5)  and  the  color  image 
corrected  by  the  single-scattering  meth¬ 
od  (Fig.  8)  reveals  that  the  oceano¬ 
graphic  features  shown  in  these  earlier 
figures  are  also  evident  in  the  third 
principal  component  image.  Thus,  the 
third  component  is,  in  effect,  an 
atmospherically  corrected  image  showing 
ocean  color  variability  patterns.  It 
should  be  noted  that  the  third  component 
will  not  always  be  the  component  to  give 
maximum  definition  of  oceanic  patterns. 
In  other  cases  where  PCA  was  applied  to 
a  small  subsection  of  an  image  so  that 
angular  variation  of  path  radiance  was 
small,  the  second  component  seemed  to 
contain  the  chlorophyll-a  variance  and, 
hence,  the  best  oceanic  patterns.  In 
general,  it  cannot  be  stated  which 
component  will  best  show  the  ocean 
patterns . 


3.3  Corrected  Imagery 

The  a,  .  coefficients  derived  from  PCA 
kj 

within  the  boxed  area  in  Figure  4  can 
now  be  used  as  a  transformation  matrix 
to  be  applied  to  the  entire  image 
according  to  Eq.  (6)-  If  the  initial 
assumptions  about  the  statistical  prop¬ 
erties  of  the  data  are  valid,  the  first 
principal  component,  which  contains 
96.6%  of  the  total  variance  in  the  sub- 
area,  can  be  assumed  to  represent  the 
path  radiance  resulting  from  aerosol 
backscatter.  For  a  check  on  this  assump¬ 
tion,  Figure  9  shows  the  first  principal 
component  and  the  original  channel  4 
image  with  the  Rayleigh  radiance  contri¬ 
bution  removed  (i.e.,  the  aerosol 


The  second  and  fourth  components  contain 
some  oceanographic  information  but  are 
predominated  by  other  influences.  The 
second  component,  with  2.8%  of  the  total 
variance,  seems  to  represent  the  angular 
dependence  of  the  atmospheric  path  radi¬ 
ance,  while  the  fourth,  with  <0.1%  of 
the  total  variance,  contains  mainly  ran¬ 
dom  noise.  For  sake  of  completeness, 
components  two  and  four  are  shown  in 
Figure  11. 

3.4  Atmospheric  Correction  Summary 

The  Grand  Banks  data  set  has  been  used 
to  show  how  a  statistical  approach  based 
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Figure  10.  The  third  principal  component  of  the  CZCS  image  shown  in  Figure  4 


on  PCA  can  be  used  to  develop  CZCS  aero¬ 
sol  correction  terms  that  are  essential¬ 
ly  equivalent  to  those  derived  from 
single- scattering  theory.  The  example 
shows  that  the  statistical  analysis  can 
be  performed  in  a  relatively  straight¬ 
forward  manner  that  does  not  have  to 
consider  such  problems  as  the  determina¬ 
tion  of  an  Angstrom  exponent,  the  valid¬ 
ity  of  ocean  blackness  at  670  nm,  sensor 
calibration,  or  many  of  the  other  prob¬ 
lems  associated  with  using  the  modeling 
approach .  Further,  the  example  shows 
that  the  statistical  approach  has 
resulted  in  an  aerosol-corrected  image 
that  apparently  shows  ocean  features.  It 
is  concluded,  therefore,  that  PCA  may  be 


an  easy  to  use  alternative  to  single¬ 
scatter  modeling  of  the  atmosphere  for 
those  applications  where  image  pattern 
interpretation,  rather  than  ocean  spec¬ 
tral  radiance  measurement,  is  the  objec¬ 
tive.  To  further  show  the  utility  of 
PCA,  four  case  studies  are  presented  in 
the  following  section  for  three  ocean 
regions . 

4.  CZCS  Case  Studies 

One  of  the  questions  posed  in  the  intro¬ 
duction  was.  Do  deep  ocean  fronts  typi¬ 
cally  have  color  signatures?'’  As  an  aid 
to  answering  this  question,  four  data 
sets  representative  of  three  different 


Figure  11. 
color  (Fig. 


The  second  and  fourth  principal  components  of  the  image  shown  in 
4)  on  page  5. 
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ocean  areas  have  been  selected  to  be 
atmospherically  corrected  by  the  PCA 
method.  These  atmospherically  corrected 
images  are  compared  with  coincident 
NOAA-IR  imagery  and  surface  data  in  an 
attempt  to  establish  the  relationship 
between  the  patterns  found  in  the  PCA 
corrected  CZCS  imagery  and  actual  ocean 
features.  The  four  test  cases  are: 

•  Case  I.  Grand  Banks,  15  May  1979; 

•  Case  II.  Gulf  of  Mexico  (summer), 

19  June  1979; 

•  Case  III.  Gulf  of  Mexico  (winter), 

6  December  1978; 

•  Case  IV.  Mediterranean,  12  October 
-1982. 

4.1  Case  I — Grand  Banks 

The  Grand  Banks  area  of  the  North  Atlan¬ 
tic  received  intensive  study  during  the 
Grand  Banks  Experiment  in  1978  and  1979 
(La  Violette,  1981).  The  Nimbus-7  CZCS 
image  chosen  for  analysis  was  acquired 
on  15  May  1979  during  the  ’’New  Look” 
portion  of  the  Grand  Banks  Experiment. 
This  case  is  accompanied  by  extensive 
surface  data. 

4.1.1  The  Surface  Data 

The  study  area  for  the  Grand  Banks 
Experiment  is  shown  in  Figure  12.  The 
bathymetric  features,  such  as  the  New¬ 
foundland  Ridge  and  the  Newfoundland 
Seamounts,  are  important  controls  of  the 
currents  in  the  area.  A  study  of  satel¬ 
lite  imagery  for  the  5-year  period, 
1975-1979,  showed  that  the  fronts  and 
eddies  associated  with  the  interaction 
of  the  Gulf  Stream,  Labrador  Current, 
and  North  Atlantic  Gyre  were  consistent¬ 
ly  positioned  in  the  area  of  the  three 
nodes  shown  in  the  simplified  schematic 
in  Figure  13  (La  Violette,  1981). 

As  might  be  expected  from  the  fact  that 
these  nodal  positions  are  correlated 
with  the  regional  bathymetry,  the  ther¬ 


mal  features  seen  in  the  satellite 
imagery  are  not  simply  surface  effects 
but  are  the  surface  manifestations  of 
dynamic  structures  that  extend  to  1000  m 
in  depth.  Figure  14  is  a  temperature 
cross-section  across  node  1  derived  from 
a  line  of  Airborne  Expendable  Bathy¬ 
thermographs  (AXBT)  recorded  on  16  May 
1979  (the  day  after  the  images  of  Fig¬ 
ures  4  and  7  were  acquired) .  This  cross 
section  clearly  shows  that  the  depth  of 
the  structure  associated  with  node  1 
extended  deeper  than  the  400-m  limit  of 
the  AXBTs . 

4.1 .2  The  CZCS  Data 

The  CZCS  data  set  for  this  case  study  is 
the  same  that  was  used  to  demonstrate 
the  atmospheric  correction  techniques  in 
Section  3.  In  that  case,  the  image  shown 
covered  an  east-west  swath  of  approxi¬ 
mately  15°  of  longitude  that  included 
nearly  the  entire  CZCS  scan  width.  How¬ 
ever,  this  discussion  covers  only  the 
eastern  half  of  the  image  (i.e.,  the 
area  located  over  nodal  features  1  and 
2).  Figure  15  is  a  color  image  formed 
from  channels  1,  2,  and  3  of  CZCS  orbit 
■2811  data..  The  scale  of  Figure  15  is 
such  that  the  spatial  resolution  is 
about  1.5  km/pixel  compared  with  the  0.8 
km/pixel  nadir  resolution  of  the  origi¬ 
nal  data.  The  Labrador  Front  nodal  posi¬ 
tions  are  marked  for  reference  purposes, 
to  show  that  the  atmospherically  uncor¬ 
rected  data  does  not  reveal  the  presence 
of  these  oceanographic  features.  Figure 
16  is  the  channel  6  thermal-IR  image 
corresponding  to  Figure  15.  Note  the 
predominance  of  the  two  nodal  features 
in  the  sea  surface  temperature  image. 

The  CZCS  image,  atmospherically  cor¬ 
rected  by  means  of  PCA  is  shown  in  Fig¬ 
ure  17.  Now,  oceanographic  patterns  are 
displayed  that  closely  resemble  those  in 
the  infrared  image  of  Figure  16,  thus 
indicating  that  the  currents  that  give 
rise  to  the  surface  thermal  patterns 
exhibit  associated  ocean  color  signa¬ 
tures  . 
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Figure  12.  Grand  Banks  Experiment  study  area.  The  square  represents  the  major 
ship  and  aircraft  survey  area  (After  La  Violette,  1982). 


4.1.3  Discussion  of  Case  I 


Although  the  general  nature  of  the  color 
and  thermal  patterns  in  Figures  16  and 
17  are  similar,  there  are  noticeable 
differences  that  are  of  oceanographic 
interest.  For  example,  the  IR  image 
shows  a  single  well-defined  front  on  the 
southwestern  side  of  node  1.  In  the 


color  image  two  fronts  are  present;  one  # 

corresponds  exactly  with  the  temperature 
front,  and  the  second  lies  parallel  to 
the  first  front  but  offset  to  the  south¬ 
west  by  40  to  50  km.  The  thermal  imagery 
does  not  indicate  a  thermal  feature  at 
the  location  of  the  second  color  front.  # 

Discussion  of  this  particular  color/tem¬ 
perature  frontal  structure  has  been 
given  by  Mueller  and  La  Violette  (1980). 
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Figure  13.  Samples  of  satellite  infrared  imagery  for  the  Grand  Banks  Experiment 
area  for  the  period  January  1975  through  October  1979.  These  imagery  show  that 
three  frontal  extrusions  were  always  present  in  the  cloud-free  data.  However, 
imagery  showing  all  three  structures  at  one  time  are  rare.  The  ones  presented  here 
are  meant  only  as  examples  of  these  features.  Newfoundland  is  the  dark  land  mass 
in  the  upper  left  corner  of  each  image.  The  line  drawing  on  the  bathymetry  chart  is 
a  composite  of  all  the  imagery  for  1978,  and  is  drawn  by  hand  on  a  common  grid.  The 
solid  line  represents  the  edge  of  the  cold  water  of  the  Labrador  Front.  The  dashed 
line  represents  the  direction  and  type  of  extrusion  away  from  the  front.  The  shaded 
area  northwest  of  the  front  represents  the  dominant  position  of  the  Labrador  Current 
as  seen  in  the  imagery.  The  smaller  shaded  areas  in  the  west  are  smaller  slope  water 
current  features  that  also  appeared  in  the  1978  imagery.  Composites  drawn  from  other 
years'  imagery  show  similar  results.  (After  La  Violette,  1982) 
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Figure  14.  Vertical  analysis  in  °C  of  data  from  aircraft  XBTs  dropped  on  16  and 
17  June  1979. 


The  oceanographic  interpretation  of  the 
ocean  color  patterns  observed  in  Figure 
17  is  conjectural.  Even  with  the  exten¬ 
sive  ground  truth  available,  one  cannot 
really  give  more  than  hypothetical 
explanations  for  what  is  observed.  How¬ 
ever,  for  the  purpose  of  this  study,  it 
is  important  to  state  that  Figure  17 
shows  that  the  interacting  current  sys¬ 
tems  in  the  Grand  Banks  area  do  have 
color  signatures  discernable  in  the  CZCS 
data . 


4.2  Case  II — Gulf  of  Mexico  (Summer) 

The  second  case  illustrates  the  effect 
of  mid-latitude  summer  atmospheric  con¬ 
ditions  (19  June  1979)  in  the  CZCS  data. 
This  example,  although  not  accompanied 
by  good  surface  data,  is  a  good  one  for 
the  purposes  of  this  study.  In  the 
image,  the  unusual  late  season  outbreak 
of  cool,  dry  continental  air  over  the 
northeastern  Gulf  of  Mexico  permits 
localized  NOAA  IR  imaging  of  sea  surface 
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Figure  15.  Case  I  Grand  Banks  CZCS  image  formed  from  channels  1,  2,  and  3 
data  uncorrected  for  atmospheric  effects.  Position  of  nodes  1  and  2  from  Figure 
1  3  are  marked  . 


Figure  16.  CZCS  channel  6  (thermal  IR)  image  corresponding  to  Figure  15. 


thermal  patterns  that  is  not  normally 
possible  in  the  Gulf  of  Mexico  in  late 
June  because  of  the  seasonal  warm,  humid 
atmosphere.  The  main  CZCS  verification 
in  this  case,  then,  will  be  a  comparison 
of  the  CZCS  color  patterns  and  the  local 
mesoscale  ocean  features  seen  in  the 
NOAA  image.  It  will  also  be  of  interest 
to  see  the  color  patterns  discernible  in 
those  areas  where  the  warm,  humid  atmos¬ 
pheric  conditions  obscure  the  ocean 
thermal  patterns. 

4.2.1  The  Surface  Data 

Although  there  is  no  specific  surface 
data  to  support  an  analysis  of  this  CZCS 


data  set,  the  NOAA-NESS  Gulf  Stream 
analysis  for  16-20  June  1979  (Fig.  18) 
does  have  the  Loop  Current  boundary 
well-defined.  According  to  NOAA-NESS 
practice,  all  available  concurrent  ship 
reports,  XBTs ,  or  AXBTs  available  from 
the  Gulf  of  Mexico  commonly  are  incorpo¬ 
rated  into  their  analysis.  So,  although 
the  accuracy  of  the  NOAA-NESS  Loop  Cur¬ 
rent  boundary  location  is  not  known,  it 
can  be  assumed  that  the  location  plotted 
on  the  chart  is  not  in  conflict  with 
available  surface  data. 

Figure  19  shows  the  TIROS-N  IR  (channel 
4)  image  from  orbit  3510  on  19  June  1979 
(the  figure  is  an  enlargement  of  a  sec- 


tion  of  Fig,  3).  Note  the  effect  of  the 
cool,  dry  air  over  the  northeastern 
Gulf.  Here  the  position  of  the  Loop 
Current  is  clearly  visible.  The  Loop 
Current,  however,  fades  from  view  to  the 
south  and  west  under  the  typical  Gulf  of 
Mexico  summer  atmosphere. 

4.2.2  The  CZCS  Data 

Figure  20  is  a  color  image  formed  from 
CZCS  channels  1,  2,  and  3  of  Nimbus-7 

orbit  3310  on  19  June  1979.  (The  image 
has  been  subsampled  to  approximately  a 
3.3  km/pixel  spacing  to  contain  the  en¬ 
tire  Gulf  of  Mexico.)  The  wide,  bright 
band, lying  in  an  approximate  northwest- 
southeast  orientation  of  the  image,  is 
apparently  caused  by  sun  reflecting  off 
the  ocean  surface.  On  this  particular 
orbit,  the  CZCS  was  tilted  to  look  be¬ 
hind  the  spacecraft  at  a  20°  tilt  angle. 
Thus,  for  this  ascending  trajectory,  the 
scanner  was  oriented  toward  the  south 
where  it  received  the  sea’s  reflection 


of  the  sun  at  its  near-zenith  summer 
time  position.  Because  of  the  presence 
of  this  sun  glint,  the  image  is  uncor- 
rectable  by  the  single-scattering  model. 
The  example  is  of  additional  interest  as 
it  demonstrates  the  advantage  of  the  PCA 
atmospheric  correction  method  for  images 
contaminated  by  sun  glint. 

The  image  subarea  chosen  for  PCA  purpos¬ 
es  is  enclosed  in  the  box  superimposed 
on  the  image  in  Figure  20.  Statistics  of 
the  data  within  the  box  are  shown  below. 

Case  II  Gulf  of  Mexico 
PCA  Box  Statistics 


mean  =  193.1,  188.9,  198.9,  238.6 


std  dev  =  3.74,  3.80, 
eigenvalues  =  0.581,  i 

0.045  0.229 

a  =  -0.517  -0.574 

3  -0.849  0.284 

0.090  -0.731 


4.70,  7.12 
.338,  0.075,  0.003 

0.341  0.910 

-0.517  0.365 

0.405  -0.181 
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Figure  18.  NOAA-NESS  Gulf  Stream  Analysis  chart  for  the  period  16-20  June  1979. 


Figure  17.  Atmospherically  corrected  version  of  Figure  15.  Red,  green,  and  blue 
£  components  of  this  image  come  from  the  second,  third,  and  second  minus  third 

principal  components  respectively. 


Table  1  shows  there  are  some  notable 
differences  between  the  data  from  this 
example  and  the  Grand  Banks  CZCS  data. 
In  this  case,  the  first  (aerosol)  prin^ 
cipal  component  contained  only  58.1%  of 
the  total  variance  within  the  PCA  box 
compared  to  Case  I  where  the  aerosol 
contribution  was  96.6%  of  the  total. 
This  variance  difference  probably 
results  from  the  combined  effects  of 
less  atmospheric  variability  and  more 
upwelled  ocean  radiance  in  the  Gulf  of 
Mexico  data.  A  similarity  between  cases 
I  and  II  is  that  the  third  component  in 
both  cases  seems  to  contain  chlorophyll 
concentration  information.  This  conclu¬ 
sion,  as  discussed  in  Section  3,  is 
based  upon  the  large  negative,  near 
zero,  and  large  positive  weights  for 
channels  1,  2,  and  3,  respectively,  that 
go  into  forming  the  third  component. 


Figure  21  shows  the  first  principal 
component  that,  under  the  assumptions  of 
the  PCA  method,  is  the  atmospheric  aero¬ 
sol  contribution  to  Figure  20.  Note  that 
the  sun  glint  band  in  the  cetiter  of  Fig¬ 
ure  20  is  present  in  the  first  principal 
component.  From  this  it  appears  that  the 
PCA  atmospheric  correction  also  removes 
some  sun  glint.  Figure  22  shows  this  is 
true.  This  figure  (the  third  principal 
component  of  Figure  20)  shows  a  complete 
absence  of  sun  glint  while  still 
displaying  what  is  believed  to  be  the 
regional  chlorophyll  distribution. 

4.2.3  Discussion  of  Case  II 

The  most  important  oceanic  features  of 
Figure  22  are  the  three  distinct  eddies, 
one  cold  (counterclockwise  rotation)  and 
two  warm  (clockwise  rotation) .  Eddies  in 


Figure  19.  TIROS-N  thermal-IR  (11.0/1)  imagery  of  the  Gulf  of  Mexico  on  19  June 
1979  at  0850Z .  This  is  the  same  image  shown  as  Figure  4,  but  at  a  different  scale. 
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the  western  Gulf  of  Mexico  are  observed 
routinely  on  IR  satellite  images  during 
the  winter  months.  However,  the  image 
shows  to  the  authors1  knowledge,  the 
first  summer  satellite  observation  of 
eddies  in  the  central  or  western  Gulf  of 
Mexico.  Note  that  the  TIROS-N  IR  image 
for  this  date  (Figure  19)  does  not  show 
these  eddies  or  the  Loop  Current  inflow 
through  the  Yucatan  straits  (also  visi¬ 
ble  in  Figure  22).  This  is  an  excellent 
example  of  the  ability  of  the  CZCS  to 
image  mesoscale  ocean  features  through  a 
warm,  humid  atmosphere  where  infrared 
images  would  not  show  the  features. 

A  color  image  formed  from  the  second, 
third,  and  fourth  principal  components 
of  the  Case  II  CZCS  image  is  shown  as 
Figure  23.  Note  that  the  Caribbean 
waters  south  of  the  Yucatan  straits  and 
the  Loop  Current  inflow  are  in  shades  of 
green,  while  the  coastal  and  shelf  wa¬ 
ters  of  the  Gulf  of  Mexico  are  in  vary¬ 
ing  shades  of  blue,  purple,  and  red.  One 
might  speculate  that  the  greenish  waters 
within  the  central  Gulf  are  of  Loop  Cur¬ 
rent  origin,  while  the  blue/purpl e/ red 
waters  are  of  Gulf  of  Mexico  origin 
(river  inflow  for  example).  The  colors 
in  Figure  23  may,  therefore,  possibly  be 
used  as  indicators  of  water  origin, 
i.e.,  tracers  for  general  circulation 
within  the  Gulf.  If  this  is  true,  Figure 
£3  would  indicate  that  most  of  the  water 
in  the  west-central  Gulf  of  Mexico  is  of 
Loop  Current  origin.  This  is  consistent 
with  predictive  models  of  the  Gulf  of 
Mexico  (Hurlburt  and  Thompson,  1980) 
that  indicate  that  eddies  shed  from  the 
Loop  Current  move  toward  the  west- 
southwest  . 

Note  that  the  two  eddies  in  the  north 
central  Gulf  are  positioned  along  the 
loop  water/ shelf  water  boundary.  Inter¬ 
pretation  of  the  color  variations  indi¬ 
cate  that  shelf  water  is  entrained  along 
the  northern  edge  of  both  the  warm 
eddies,  and  Loop  Current  water  is  en¬ 
trained  at  the  southern  edges.  These 
eddies,  in  effect,  appear  to  be  mixing 
areas  for  the  two  contrasting  water 


types.  If  this  is  true,  some  interesting 
speculation  may  be  made.  Hawkins  (1983) 
has  observed  very  rapid  surface  cooling 
of  eddies  separated  from  the  Loop  Cur¬ 
rent.  From  considerations  of  heat  con¬ 
tent  and  possible  heat  fluxes  through 
the  sea/air  boundary,  the  observed  cool¬ 
ing  rates  seem  to  be  unrealistically 
large.  Entrainment  and  subsequent  mixing 
of  substantial  amounts  of  shelf  water 
(as  the  study  of  Fig.  23  suggests)  may 
explain  the  rapid  decay  of  eddy  surface 
thermal  signatures  in  this  area.  Thus, 
the  CZCS  may  be  a  tool  to  disclose  mix¬ 
ing  processes  in  the  ocean. 

4.3  Case  III — Gulf  of  Mexico  (Winter) 

The  Remote  Sensing  Branch  of  NORDA  con- 
ducted  a  cruise  in  the  Gulf  of  Mexico  in 
December  of  1978.  CZCS  data  from  a  por¬ 
tion  of  that  cruise  period  (2  December 
1978),  has  been  chosen  as  a  Gulf  of 
Mexico  wintertime  case  study.  These  data 
will  be  compared  to  XBT,  chlorophyll, 
and  other  water  color  related  measure¬ 
ments  collected  during  the  NORDA  cruise. 


4.3.1  The  Suface  Data 

The  IR  imagery  used  for  comparison  with 
the  CZCS  visible  imagery  for  this  case 
is  derived  from  the  CZCS  IR  channel  of 
the  same  CZCS  data  set.  Although  the 
CZCS  thermal  images  are  crude  by  compar¬ 
ison  with  the  NOAA  imagery,  clouds  pres¬ 
ent  during  the  orbits  of  N0AA-5  and 
TIROS-N  precluded  their  use,  leaving  the 
CZCS  orbit  the  only  cloud-free  imagery. 

The  line  segment  AB  in  Figure  25  marks 
an  XBT  line  made  from  the  USNS  LYNCH  on 
6  December  1978,  (4  days  after  the  CZCS 
pass).  Figure  24  is  a  vertical  tempera¬ 
ture  cross  section  plotted  from  those 
XBTs .  By  convention  the  Loop  Current 
boundary  is  considered  to  be  where  the 
22°C  isotherm  drops  below  a  depth  of  100 
m  (Leipper,  1970),  and  from  this  the 
line  segment  AB  in  Figure  24  straddles 
the  Loop  Current  boundary.  Therefore, 
the  warm  (dark)  feature  seen  in  the 


Figure  20.  CZCS  image  formed  from  bands  1,  2,  and  3  of  orbit  3310  over  the 
Gulf  of  Mexico  on  19  June  at  1729  Z.  The  red  box  marks  the  area  included  in 
the  PCA  analysis. 


Figures  21  and  22.  First  principal  component  (top)  of  the  Case  II  CZCS  image  shown 
in  Figure  20.  The  third  principal  component  (bottom,  primarily  showing  chlorophyll 
concentration)  of  the  CZCS  image  shown  in  Figure  20. 
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Figure  23.  Atmospherically  corrected 
image  formed  from  the  second,  third. 


(RCA  method)  version  of  Figure  20. 
and  fourth  principal  components. 
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Figure  24.  XBT  data  from  along  line  AB  in  Figure  25  shown  as  contours  of  tempera¬ 
ture  as  a  function  of  position  and  depth. 


#  south-eastern  quadrant  of  the  Gulf  is 
indeed  the  Loop  Current.  Furthermore, 
the  XBTs  show  that  the  warm  feature  in 
the  IR  imagery  is  not  simply  a  surface 
effect,  but  is  associated  with  oceanic 
structure  that  extends  at  least  to  the 

#  depth  range  of  the  XBTs. 

Chlorophyll-a  concentration  and  total 
suspended  solids  measurements  were  made 
on  both  sides  of  the  Loop  Current  bound¬ 
ary  during  the  USNS  LYNCH  cruise.  These 

#  data  show  that  surface  chlorophyll-a 
concentrations  within  the  Loop  Current 
on  5  December  at  point  A  of  the  line 
segment  AB  in  Figure  25  was  0.05  mg/1  at 
the  surface  and  0.30  mg/1  at  20  m.  Out¬ 
side  the  Loop  Current  at  point  C  in  Fig- 

®  ure  25,  the  chlorophyll— a  concentration 

at  depths  of  2  and  20  m  was  0.19  and 
0.14  mg/1  when  measured  on  12  December. 


Although  these  two  measurements  were  7 
days  apart  and  neither  coincided  with 
the  CZCS  data  on  2  December,  these  data 
show  that  a  large  difference  in  the  ver¬ 
tical  distribution  of  chlorophyll-a  did 
exist  between  Loop  Current  and  Gulf 
waters.  Although  the  values  on  2  Decem¬ 
ber  may  have  been  different,  these 
measurements  indicate  that  a  strong 
chlorophyll-a  change  did  take  place 
across  the  Loop  Current  boundary.  Total 
suspended  solids  at  these  stations 
(average  of  surface  and  20-m  depth  val¬ 
ues)  were  490  ug/1  at  point  A  and  195 
ug/1  at  point  C  indicating  that  solids 
other  than  chlorophyll-a  were  also  in 
the  water  in  varying  concentrations. 
These  measurements  would  suggest  that 
one  could  expect  to  see  a  color  differ¬ 
ence  between  the  resident  Gulf  of  Mexico 
and  Loop  Current  waters. 
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Figure  25.  CZCS  channel  6  (thermal  IR)  from  orbit  545  over  the  Gulf  of  Mexico  on  2  December  1978  at  1715Z. 
The  line  AB  marks  the  location  of  an  XBT  transect  by  the  USNS  LYNCH  on  6  December  1978. 


Figure  26.  CZCS  image  formed  from  bands  1,  2,  and  3  of  orbit  545  over  the  Gulf  of  Mexico 
on  2  December  1978  at  1715  Z.  The  red  box  marks  the  area  included  in  the  PCA  analysis. 


4.3.2  The  CZCS  Data 

Figure  26  is  a  color  image  formed  from 
channels  1,  2,  and  3  of  CZCS  orbit  545 
on  2  December  1978.  This  image  has  been 
transformed  to  Mercator  projection  with 
the  same  3.3  km/pixel  scale  as  Case  II. 
The  image  subarea  chosen  for  PCA  analy¬ 
sis  is  indicated  by  the  red  box  in  the 
figure.  Data  statistics  within  the  box 
are  as  follows : 

Case  III  Gulf  of  Mexico 
PCA  Box  Statistics 

mean  =  168.6,  133.6,  131.0,  116.7 


std 

dev  =  1.52, 

1.73, 

2.31 

y  4 « 

,47 

eigenvalues  =  0 

.893,  0 

.069 

,  0. 

,023 

,  0.1 

0.180 

0.293 

0. 

410 

0. 

844 

=  -0.723 

-0.371 

-0. 

358 

0. 

458 

kj 

-0.616 

0.167 

0. 

718 

-0. 

275 

0.252 

-0.864 

0. 

432 

0. 

036 

As 

in  the  previous  cases 

the 

akj 

mat 

shows  a  component  (the  third  one  in  this 
case)  that  seems  to  capture  chlorophyll- 
a  variability  with  its  -0.616,  -0.167, 
and  0.718  weights  in  channels  1,  2,  and 
3  respectively.  Figure  27  shows  the 
first  principal  component  that,  under 
the  assumptions  of  the  PCA  method,  is 
the  atmospheric  aerosol  contribution  to 
Figure  26. 

4.3.3  Discussion  of  Case  III 


An  atmospherically  corrected  Case  III 
CZCS  color  image  (Fig.  28)  formed  from 
principal  components  2,  3,  and  4  shows 
the  Loop  Current  boundary  to  be  in  the 
same  location  as  that  shown  by  XBT  data 
in  Figure  24.  The  position  and  shape  of 
the  Loop  Current  pattern  in  the  color 
image  corresponds  with  the  position  and 
shape  of  the  warm  feature  in  the  thermal 
image  (Fig.  24)  indicating  that  the 
color  scanner  is  indeed  imaging  the  Loop 
Current.  The  Loop  Current  is  the  only 
major  deep-ocean  feature  apparent  in 
Figure  28.  This  is  consistent  with  the 
finding  of  the  USNS  LYNCH  data  that 


showed  that  no  eddies  or  major  oceano¬ 
graphic  features,  other  than  the  Loop 
Current,  were  found  in  the  eastern  Gulf 
at  the  time. 

The  main  conclusion  to  be  drawn  from 
Case  III  is  that  the  Loop  Current 
appears  to  have  a  color  signature  in  the 
wintertime  as  well  as  in  the  Case  II 
summertime  image.  This  would  tend  to 
indicate  a  persistence  in  the  color 
uniqueness  of  the  Loop  Current  through¬ 
out  the  year.  This  is  important  from  a 
Navy  point  of  view  since  if  the  Loop 
Current,  or  other  similar  mesoscale  fea¬ 
tures,  possessed  an  observable  color 
signature  only  during  short-lived  tran¬ 
sient  biological  events,  the  utility  of 
ocean  color  for  Navy  operations  would  be 
diminished . 

4.4  Case  IV:  Mediterranean 

NORDA  scientists  participated  in  an 
international  oceanographic  experiment , 
called  6Donde  Va? ,  conducted  in  the 
western  Mediterranean  in  October  1982 
(La  Violette,  et  al.  1982;  Kinder, 
1983).  The  study  focused  on  the  Alboran 
Gyre,  a  large  eddy  found  in  the  western 
Alboran  Sea.  The  Alboran  Sea  is  an 
oceanographically  dynamic  area  where 
Mediterranean,  Atlantic,  locally  up- 
welled,  and  coastal  waters  constantly 
interact.  This  area,  therefore,  tends  to 
produce  a  wide  variety  of  ocean  color 
signatures  that  can  be  easily  sensed  by 
the  CZCS. 

4.4.1  The  Surface  Data 

During  the  period  6-18  October  1982, 
ship,  aircraft,  and  satellite  data  were 
collected  as  part  of  the  iDonde  Va? 
experiment  in  the  Alboran  Sea  (La 
Violette,  1983).  The  ship  data  included 
CTD,  XBT,  and  chlorophyll  measurements, 
whereas  the  aircraft  data  included  sea 
surface  temperature,  sonobuoy  drift,  and 
wind  measurements.  The  satellite  data, 
acquired  every  12  hours,  included  infra¬ 
red  and  visible  measurements  that  were 
registered  and  atmospherically  corrected 
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Figure  27.  The  first  principal  component  (aerosol)  of  the  Case  III  CZCS  image. 


Figure  28.  Color  image  formed  from  the  second,  third,  and  fourth  principal  components  of  the 
Case  III  CZCS  image. 


0349  GMT  11  OCTOBER  1982 


1513  GMT  11  OCTOBER  1982 


0337  GMT  12  OCTOBER  1982 


1501  GMT  12  OCTOBER  1982 


1449  GMT  13  OCTOBER  1982 


Figure  29.  NOAA-7  infrared  imagery  for  the  early  morning  of  11  October 
through  the  afternoon  of  13  October  1982. 
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to  form  sea  surface  radiation  images  for 
the  study  period.  Analysis  of  these 
measurements  shows  that  incoming  Atlan¬ 
tic  water,  as  it  left  the  Strait  of 
Gibraltar,  moved  at  the  surface  in  a 
narrow  (approximately  16  km  wide)  band 
of  approximately  17-18°C  water  that 
formed  the  initial  quarter  of  the  cir¬ 
cumference  of  the  gyre.  As  the  water 
continued  the  incoming  water  was  modi¬ 
fied;  however,  it  still  retained  its 
identity.  Associated  with  the  inflowing 
band  of  water  were  cold  water  features 
that  hung  southward  toward  the  center  of 
the  gyre.  These  also  moved  around  the 
gyre,  changing  shape  as  they  moved. 


The  five-image  sequence  of  Figure  29 
shows  the  movement  of  several  of  these 
mesoscale  cold-water  features  around  the 
gyre.  Notice  the  feature  developing  on 
12  October  in  the  region  east  of  Gibral¬ 
tar.  An  enlargement  of  the  IR  image  for 
the  afternoon  of  12  October  (1501  hrs 
GMT)  is  presented  in  Figure  30  together 
with  the  principal  component  for  chloro¬ 
phyll  from  the  Nimbus-7  CZCS  visible 
image  of  the  same  area  for  1155  hrs  GMT 
of  the  same  day.  Ship  data  taken  for  the 
same  period  show  the  pattern  of  distri¬ 
bution  of  the  chlorophyll  to  be  similar 
to  that  shown  by  the  satellite  imagery. 
Examination  of  the  NOAA  imagery  shows 
that  the  predominant  feature  moved 
approximately  17  km  in  the  11  hours  from 
the  time  of  the  morning  NOAA  IR  image  to 
the  time  of  the  afternoon  NOAA  IR  image, 
or  at  an  average  speed  of  approximately 
0.4  m/sec.  The  CZCS  image  indicates  a 
movement  of  approximately  9  km  in  the 
time  between  the  Nimbus— 7  noon  pass  and 
the  NOAA  afternoon  pass  (i.e.,  if  one 
assumes  that  the  sharpest  gradient 
change  in  the  CZCS  image  is  directly 
associated  with  the  sharpest  thermal 
gradient  in  the  NOAA  imagery) .  This 
equates  to  a  movement  of  0.8  m/sec  or 
twice  the  rate  of  movement  for  the  11- 
hour  period  between  the  two  NOAA  images. 
The  figure  gives  a  graphic  indication  of 
how  much  movement  can  take  place  in  4 
hours.  The  possible  reason  for  the  vari¬ 
ations  in  speed  is  the  phase  of  the 


diurnal  tide.  Comparison  of  the  shifts 
in  speed  with  the  tide  tables  from  Ceuta 
show  the  two  are  in  phase. 

4.4.2  The  CZCS  Data 

Eleven  essentially  cloud-free  CZCS 
images  were  collected  for  the  iDonde  Va? 
investigations  during  the  first  3  weeks 
of  October  1982.  The  12  October  image 
chosen  for  analysis  here  is  the  poorest 
image  of  the  dDonde  Va?  data  set  in 
terms  of  atmospheric  contamination.  Many 
of  the  images  during  the  experiment  were 
useable  for  descriptive  analysis  without 
any  atmospheric  correction,  but  the  12 
October  image  (Fig.  31)  is  so  severely 
contaminated  by  atmospheric  aerosols 
that  ocean  color  features  are  completely 
masked.  This  image  is  used  as  the  fourth 
case  study  in  order  to  evaluate  the  PCA 
atmospheric  correction  method  on  an 
extreme  case. 

The  PCA  analysis  was  performed  on  the 
area  shown  by  the  red  box  in  Figure  31. 
Statistics  for  the  PCA  box  are  as 
follows : 


Case  IV  Mediterranean 
PCA  Box  Statistics 


mean  =  107.8,  111.1,  117.9,  129.4 
std  dev  «  3.47,  6.13,  8.15,  15.62 
eigenvalues  =  .986,  .008,  .004,  .001 
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0.428 

-0.594 

-0.198 

0.651 


0.827 

0.557 
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Figure  32  is  the  first  principal  compo¬ 
nent  (atmospheric  aerosol)  of  the  Case 
IV  CZCS  image.  Figure  33  is  a  false 
color  image  produced  by  principal  compo¬ 
nents  2,  3,  and  4  that  contribute  red, 
green,  and  blue  to  the  color  image.  Fig¬ 
ure  33  represents  the  atmospherically- 
corrected  Case  IV  CZCS  image.  The  anti- 
cyclonic  circulation  of  the  Alboran  gyre 
is  seen  in  Figure  33  in  shades  of  gold. 
Patterns  created  by  turbid  coastal  wa¬ 
ters  are  seen  in  shades  of  purple  along 
the  African  coast. 
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Figure  30.  The  short  term 
movement  of  feature  shown  by 
Nimbus  CZCS  and  NOAA-7 
AVHRR-IR  imagery.  The 
movement  took  place  between 
1155  (Nimbus)  and  1501  (NOAA) 
hrs  GMT  on  12  October. 
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Figure  31.  CZCS  image  formed  from  bands  1,  2,  and  3  of  orbit  20029  over  the  Mediterranean 
Sea  on  12  October  1  982  at  1126  Z.  The  yellow  box  marks  the  area  included  in  the  PCA  analysis 


*  * 


Figure  32.  First  principal  component  of  the  Case  IV  CZCS  image  shown  in  Figure 
30.  Under  the  PCA  method  of  atmospheric  correction,  this  is  taken  to  be  the 
aerosol  contribution  to  the  image. 


4.4.3  Discussion  of  Case  IV 

The  PCA  method  seems  to  have  done  a 
reasonable  job  of  extracting  water  color 
features  through  this  very  hazy  atmos¬ 
phere.  Comparison  of  Figure  33  with  the 
NOAA-7  IR  image  for  the  same  day  (Fig. 

34)  shows  that  most  of  the  ocean  fea¬ 
tures  apparent  in  the  IR  image  are  also 
seen  in  the  atmospherically  corrected 
CZCS  image.  However,  some  uncorrected 
artifacts  of  the  atmosphere  do  remain  in 
Figure  33. 

The  difference  in  color  shades  for  the 
Alboran  gyre  and  the  coastal  features  in 
Figure  33  is  of  interest  in  this  case. 
It  appears  that  the  PCA  transf ormation 
has  separated  the  chlorophyll  (Alboran 
gyre)  and  inorganic  sediment  (coastal 
features)  contributions  to  the  water 
color  into  different  principal  compo¬ 
nents.  Looking  at  the  individual  princi¬ 
pal  components  that  make  up  the  image  in 
Figure  33,  we  find  that  coastal  patterns 
appear  in  the  second  component  (Fig. 

35) ,  which  does  not  show  the  gyre  at 
all.  The  gyre  appears  in  the  third  com¬ 
ponent  (Fig.  36)  that  exhibits  only  a 
very  faint  trace  of  the  coastal  fea¬ 
tures.  The  authors  have  processed  a  Gulf 
of  Mexico  data  set  (not  included  as  a 
case  study  here)  where  a  similar  separa¬ 
tion  of  chlorophyll  and  sediment  seems 
to  have  occurred. 

5.  Universal  Quick-Look  CZCS  Atmospheric 
Correction  Algorithms 

In  each  of  the  case  studies,  we  employed 
Principal  Components  Analysis  to  derive 
a  transformation  matrix  that  when 
applied  to  the  original  CZCS  spectral 
bands  produced  a  new  multiband ed  image 
with  the  atmospheric  radiance  largely 
confined  to  a  single  band;  the  other 
bands  are  predominantly  composed  of 
ocean  radiance.  The  result  of  using  this 
transformation  matrix,  therefore,  is  a 
form  of  atmospheric  correction.  The 
question  is  now  asked,  "How  similar  are 
these  transformation  matrices  from  case 


to  case?"  If  there  is  a  strong  similari¬ 
ty,  then  a  universal  set  of  matrix  coef¬ 
ficients  could  be  defined  that  could  be 
used  for  atmospheric  correction  on  most 
images  without  requiring  the  PCA  proce¬ 
dural  step. 

In  the  discussion  of  the  previous 
sections,  some  similarity  has  been 
noticeable.  All  cases  resulted  in  a 
chlorophyll-like  principal  component 
where  large  negative,  near  zero,  and 
large  positive  coefficients  for  channels 
1,  2,  and  3,  respectively,  seemed  to 

reflect  the  chlorophyll  "hinge  point" 
behavior.  The  most  striking  similarity 
observed  is  between  case  studies  I  and 
III,  Grand  Banks  and  Gulf  of  Mexico  win¬ 
ter,  where  the  "chlorophyll  component" 
weights,  shown  in  the  table  below,  are 
very  similar. 


Chlorophyll  Component  Coefficients 


CHI 

CH2 

CH3 

CH4 

Case 

I  -0.529 

0.199 

0.743 

-0.356 

Case 

III  -0.616 

0.167 

0.718 

-0.275 

The  PCA  boxes  for  these  two  cases  were 
offshore,  well  away  from  any  influence 
of  coastal  or  shelf-type  waters.  This 
type  of  location  may  explain  their  simi¬ 
larity.  If  so,  then  an  average  of  these 
two  sets  of  weights  might  be  a  good 
first  approximation  to  a  universal 
atmospheric  correction  algorithm  that 
would  result  in  delineation  of  chloro¬ 
phyll  patterns  for  deep  ocean  waters. 
Such  an  algorithm  would  be 

deep  ocean  chlorophyll  =  -  0.572*Chl  + 
0 . 183*Ch2  +  0 . 730*Ch3  -  0.315  *Ch4  (7) 

where  Chi,  Ch2,  Ch3 ,  and  Ch4  are  the  raw 
digital  counts  in  each  of  the  first  4 
CZCS  spectral  channels. 

For  case  studies  II  and  IV,  Gulf  of 
Mexico  summer  and  Mediterranean,  the 
boxes  selected  for  PCA  were  in  shelf 
waters.  These  two  cases  both  resulted  in 
"chlorophyll  components"  that  were  also 


Figure  33.  Atmospherically  corrected  (PCA  method)  version  of  Figure  31  CZCS  image  formed 
from  the  second,  third,  and  fourth  principal  components. 
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Figure  34.  NOAA-7  thermal  IR  (11.0  ju)  imagery  of  the  Alboran  Sea  on  12  October 
1982  at  1505Z. 
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Figure  35.  The  second  principal  component  for  the  Case  IV  CZCS  image.  Image 
shows  primarily  the  coastal  features  along  the  African  coast.  Evidence  of  the 
Alboran  gyre  is  absent. 


•  • 
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Figure  36.  The  third  principal  component  for  the  Case  IV  CZCS  image.  Image 
clearly  shows  the  position  of  the  Alboran  gyre. 


similar  to  each  other,  but  significantly 
different  from  the  Case  I  and  III 
results  from  PCA  in  open  ocean  waters. 
The  table  below  gives  the  "chlorophyll 
component"  weights  for  cases  II  and  IV. 


Chlorophyll  Component  Coefficients 


Chl 

Ch2 

Ch3 

Ch4 

Case 

Case 

II 

IV 

-0.849 

-0.969 

0.284 

0.137 

0.405 

0.198 

-0.181 

0.046 

These  two  sets  of  numbers  are  similar, 
but  not  as  closely  matched  as  the  Case  I 
and  III  deep  water  values.  More  varia¬ 
bility  in  shelf  areas  is  probably  to  be 
expected  because  of  the  greater  varia¬ 
bility  of  water  types. 

Here  again,  we  propose  that  the  average 
of  these  two  sets  of  coefficients  might 
be  a  first  guess  for  a  usable,  univer¬ 
sal,  quick-look  algorithm  for  imaging 
chlorophyll  distribution  in  shelf 
waters.  The  algorithm  would  be 

shelf  chlorophyll  *  -  0.909*Chl  + 

0 .210*Ch2  +  0 . 301*Ch3  -  0.067*Ch4  (8) 

How  well  would  one  expect  a  universal 
algorithm  to  perform  on  a  CZCS  image 
that  is  not  one  of  the  case  studies 
presented  here?  The  spectral  variability 
of  aerosol  phase  functions,  CZCS  radio- 
metric  calibration  uncertainties,  and 
many  of  the  same  factors  that  make  the 
single-scattering  method  difficult  to 
use  also  introduce  variability  that 
would  tend  to  degrade  the  performance  of 
any  universal  algorithms.  However,  it  is 
believed  that  the  degradation  resulting 
from  these  factors  will  not  be  suffi¬ 
ciently  large  to  prevent  their  use  for 
many  Navy  operational  or  descriptive 
applications . 

To  illustrate  this  point,  the  universal 
algorithm’s  coefficients  proposed  here 
have  been  applied  to  a  fifth  CZCS  image 
of  the  Gulf  of  Tehuantepec  off  the  west 
coast  of  Mexico.  The  original  CZCS  image 
recorded  on  30  November  1981  is  shown  in 


Figure  37.  The  results  of  the  deep  ocean 
and  shelf  water  chlorophyll  algorithms 
are  shown  as  Figures  38  and  39,  respec¬ 
tively.  These  figures  do  show  a  signifi¬ 
cant  reduction  in  atmospheric  haze,  thus 
permitting  improved  display  of  the  ocean 
mesoscale  color  patterns. 

6.  Summary  and  Conclusions 

Previous  methods  that  perform  atmospher¬ 
ic  corrections  on  CZCS  data  based  on 
single-scattering  theory  are  shown  in 
this  report  to  be  difficult  to  use.  This 
arises  from  uncertainties  concerning  the 
spectral  dependency  of  the  aerosol  phase 
function,  instrument  calibration,  as 
well  as  the  validity  of  the  single¬ 
scattering  approximation  for  optically 
thick  atmospheres.  A  statistical 
approach  based  on  Principal  Components 
Analysis  that  avoids  many  of  these  prob¬ 
lems  is  presented. 

The  statistical  method,  unlike  the  meth¬ 
ods  based  on  single-scattering  theory, 
does  not  result  in  quantitative  esti¬ 
mates  of  spectral  radiance.  The  PCA  cor¬ 
rection  method  is,  therefore,  not  pro¬ 
posed  as  a  replacement  for  the  single¬ 
scattering  methods  in  those  applications 
where  spectral  radiance  measurements  are 
needed.  The  new  method  is  useful  for 
those  cases  where  a  descriptive  inter¬ 
pretation  of  the  corrected  imagery  is 
required . 

The  four  case  studies  analyzed  seem  to 
indicate  that  mesoscale  features  in  the 
deep  ocean  are  frequently  discernable  in 
atmospherically  corrected  CZCS  imagery. 
The  case  study  of  the  Alboran  Sea  showed 
the  features  could  be  used  as  part  of 
the  analysis  of  the  movement  of  meso¬ 
scale  features.  Thus,  the  use  of  this 
type  of  data  is  promising  for  Naval 
applications  where  mesoscale  structure 
must  be  identified  and  monitored.  The 
corrected  imagery  has  been  shown  to  have 
other  uses  in  these  case  studies.  The 
summertime  Gulf  of  Mexico  case  study  has 
demonstrated  the  imaging  of  mesoscale 
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Figure  37.  Color  image  formed  from  CZCS  bands  1,  2,  and  3  on  orbit  1  5665  at  1750  Z 
on  30  November  1981.  Area  shown  is  the  Gulf  of  Tehuantepec  on  the  Pacific  Coast  of 

Mexico. 
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Figure  38.  Image  formed  from  application  of  the  universal  deep-water  algorithm 
(Eq.  7)  to  the  image  of  Figure  37. 
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Figure  39.  Image  formed  from  application  of  the  universal  shelf-water  algorithm 
(Eq.  8)  to  the  image  of  Figure  37. 
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features  through  a  warm,  humid  atmos¬ 
phere  that  was  opaque  to  IR  sensors.  The 
ability  of  the  Nimbus-7  CZCS  to  see 
through  these  atmospheres  where  IR  re¬ 
mote  sensing  of  the  ocean  is  ineffective 
is  potentially  important  to  the  Navy. 

Some  consistency  between  the  case 
studies  indicates  that  universal  algo¬ 
rithms  for  atmospheric  correction  based 
on  a  simple  weighted  sum  of  channels  may 
be  possible  without  going  through  the 
PGA  procedure  for  every  case.  While  not 
optimal,  universal  weighting  coeffi¬ 
cients  for  atmospheric  correction  by 
summing  CZCS  channels  may  be  adequate 
for  quick“look  for  Navy  operational 
purposes.  Coefficients  for  universal 
algorithms  based  on  these  four  case 
studies  have  been  proposed. 

Other  useful  techniques  have  resulted 
from  the  study.  For  example,  the  aerosol 
path  radiance  has  been  associated  with 
the  first  principal  component  of  CZCS 
images  when  the  PCA  is  performed.  The 
first  principal  component  that  results 
is  discarded  to  achieve  an  atmospheric 
correction.  However,  there  are  other 
Navy  applications  that  might  find  the 
integrated  aerosol  scattering  over  the 
ocean  as  a  desired  parameter.  In  these 
cases,  the  first  principal  component  of 
the  PCA  method  may  be  useful  in  that, 
rather  than  being  discarded  as  an 
unwanted  contribution  to  the  image,  it 
would  actually  be  the  final  product  of 
the  process  indicating  atmospheric  aero¬ 
sol  concentrations. 
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